我正在尝试编写一个函数,在二维数组中创建一个5x5像素的规则网格。我希望numpy.arange和numpy.repeat的一些组合可以做到这一点,但到目前为止我还没有任何运气,因为numpy.repeat只会沿着同一行重复。这是一个例子:假设我想在形状为(20,15)的二维数组中放置一个5x5的网格。它应该看起来像:array([[0,0,0,0,0,1,1,1,1,1,2,2,2,2,2],[0,0,0,0,0,1,1,1,1,1,2,2,2,2,2],[0,0,0,0,0,1,1,1,1,1,2,2,2,2,2],[0,0,0,0,0,1,1,1,1,1,2,2,2,2,2],[
我刚刚使用F2PY将Fortran90子例程包装到python。这里的微妙之处在于Fortran子例程aslo将python回调函数作为其参数之一:SUBROUTINEf90foo(pyfunc,a)real(kind=8),intent(in)::a!f2pyintent(callback)pyfuncexternalpyfunc!f2pyreal*8y,x!f2pyy=pyfunc(x)!***debugbegins***print*,'StartLoop'doi=1,1000p=pyfunc(a)enddototal=etime(elapsed)print*,'End:total